Constraints on Yukawa- Type Deviations from Newtonian Gravity at 20 Microns 
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Recent theories of physics beyond the standard model have predicted deviations from Newtonian 
gravity at short distances. In order to test these theories, we have a built an apparatus that can 
measure attonewton-scale forces between gold masses separated by distances on the order of 25 /an. 
A micromachined silicon cantilever was used as the force sensor, and its displacement was measured 
with a fiber interferometer. We have used our measurements to set bounds on the magnitude a and 
length scale A of Yukawa-type deviations from Newtonian gravity; our results presented here yield 
the best experimental limit in the range of A = 6-20 fim. 

PACS numbers: 04.80.Cc 



I. INTRODUCTION 



A. Motivation 



The theories of Newton and Einstein are powerful and 
accurate in describing the observed natural world. In 
spite of this, gravity still presents several theoretical chal- 
lenges, including the gauge hierarchy problem, the cos- 
mological constant problem, and the lack of a quantum 
description of gravity. Many theories of physics beyond 
the standard model, in particular those theories that at- 
tempt to unify the standard model with gravity, predict 
the existence of extra dimensions, exotic particles, or new 
forces that could cause a mass coupling in addition to the 
Newtonian gravitational potential at short distances. Of 
particular interest to us have been some recent theories 
that predict new forces in a range measurable by table- 
top experiments 0,|2>lil>EllSIE0>ll|- In many cases, 
the modification to Newtonian gravity is predicted to be 
a Yukawa-type potential, arising either from coupling to 
massive particles or the compactification of extra dimen- 
sions. With this addition, the gravitational potential be- 
tween two masses m\ and vni separated by distance r is 
predicted to be of the form: 



V(r) = -G 



m\m2 



1 + ae 



-(r/A) 



(1) 



Here, G is Newton's constant, a is the strength of the new 
potential as compared to the Newtonian gravitational po- 
tential, and A is its range. 

With these ideas in mind, we have constructed a device 
to measure attonewton-scale forces between masses se 
arated by distances on the order of 25 /jm 
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FIG. 1: [color online] Experimental results (solid lines) and 
theoretical predictions shown in a- A space. The area to the 
upper right of the experimental lines (from Refs. 0, ll4L ITM 
ITq . l3l| and our latest results from this paper) shows where 
Yukawa-type deviations from Newtonian gravity have been 
excluded. The line labeled "Stanford 1" gives results from 
0; the line labeled "Stanford 2" shows the results described 
in this paper. Dashed lines show roughly the predictions for 
the dilaton 8] and the first Kaluza-Klein mode of two simply- 
compactified extra dimensions as described in Ref. 2]. Shaded 
regions to the left of the experimental lines show predictions 
for moduli and gauge bosons from Ref. [32^ |. 



Our experiment was designed to measure Yukawa- 
type forces in the range of A = 5-50 fim with as small a 
value of I a I as possible. This range of parameters is rel- 
evant to recent theoretical predictions 0, 0, Q and is 
complementary to recent experimental attempts to test 
these new ideas 0, Hj, H3 ■ Our results are reported 
in terms of an upper bound on a (A) of a Yukawa poten- 
tial that could be consistent with our data; we present 
results at the 95% confidence level. Unless otherwise 
stated, all measurements and techniques described ap- 
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FIG. 2: [color online] Schematic (not to scale) showing the 
drive mass below the cantilever that bears the gold test mass. 
A gold-coated silicon nitride shield membrane separates the 
masses. 

ply to the Cooldown 04 experiment, from which the new 
a(X) bound was derived. These latest results yield close 
to an order of magnitude improvement over our previous 
bounds on the Yukawa potential at distances on the or- 
der of 20 /im 0. Fig. ^summarizes our results together 
with some of the theoretical predictions and other recent 
experimental results. 

B. Overview 

To measure very short-range forces, we used masses of 
size comparable to the distance between them. A gold 
prism was attached to the end of a single-crystal silicon 
cantilever. A drive mass, comprising a gold meander pat- 
tern embedded in a silicon substrate to create an alter- 
nating pattern of gold and silicon bars, was mounted on 
a piezoelectric bimorph at a distance from the cantilever- 
mounted test mass. The face-to-face vertical separation 
between the masses was ~ 25 /im, limited in part by the 
presence of a stiff, metallized silicon nitride shield mem- 
brane between the masses. 

As diagrammed in Fig. [3 the drive mass was oscillated 
along the y-direction underneath the test mass, main- 
taining the vertical (z-) separation between the masses. 
Due to the differing mass densities of the gold and sili- 
con bars, an alternating gravitational field was created at 
the test mass location. The drive mass was oscillated at a 
subharmonic of the resonant frequency of the cantilever. 
Due to the geometry of the drive mass and the ampli- 
tude of oscillation, any gravitational coupling between 
the masses would thus create a force on the cantilever 
at harmonics of this drive frequency, including the can- 
tilever's resonant frequency. The motion of the cantilever 
on resonance was measured using a fiber interferometer; 
from this measurement, the force between the masses was 
deduced. 

Thermal noise provided a limit on the measurement 
of cantilever motion. To minimize this limit, cantilevers 



were fabricated to have small spring constants and high 
quality factors in vacuum. The thermal noise limit in 
this experiment was approximately 2.5 x 10~ 16 N/vTfz, 
at cryogenic temperatures in vacuum. 

In order to be able to test and characterize the appa- 
ratus by measuring a known force much larger than any 
expected gravitational force, a magnetic analog to the 
gravitational experiment was built into the apparatus. 
Magnetic test masses were fabricated with a nickel layer 
on them. An electrical current passed through the drive 
mass meander created a spatially-varying magnetic field 
that would couple to the magnetic moment of the test 
mass and drive the cantilever. 

For both magnetic and gravitational tests, the force 
was measured as a function of the equilibrium y-position 
between the drive mass and the test mass. Any cou- 
pling between masses would show a distinct periodicity 
in the measured force as a function of this y-position. 
By comparing our measurements to predictions from fi- 
nite element analysis (FEA), a bound on Yukawa- type 
deviations was derived. 

The separation of the signal frequency from the drive 
frequency, the shield between the masses, and the use of 
non-magnetic test masses for the gravitational measure- 
ment reduced or eliminated many possible sources of non- 
gravitational background. The geometry of the appara- 
tus provided an important degree of freedom that permit- 
ted a discrimination of true coupling between the masses 
from certain electrical or mechanical backgrounds. 

C. Organization of the Paper 

This paper primarily describes the apparatus and 
results from the experiment labeled Cooldown 04. 
Cooldown 01 is described in Refs. 0,0, Cooldown 02 
is described in Sec. lVll of this paper, and Cooldown 03 is 
described elsewhere|lll Il3| . 

In the second section of this paper, the apparatus is 
described in detail. The following sections describe our 
finite element analysis, the experimental method, and the 
averaging used to convert the raw data to a force mea- 
surement. The sixth section of this paper describes the 
magnetic experiment used to test the apparatus and the 
measurement of thermal noise. In the final sections, our 
experimental results, fitting techniques, error analysis, 
and the resulting bound on a(A) are presented, followed 
by conclusions and a discussion of future prospects for 
this experiment. 

II. APPARATUS 

A. Cantilever 

The cantilevers used in this experiment were fabricated 
from single-crystal silicon using standard micromachin- 
ing techniques |10| . The cantilevers were fabricated from 
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(100) Si, oriented in the (110) direction; they were 50 /im 
wide, 250 /im long, and 0.33 /im thick, yielding an ex- 
pected spring constant k = 0.005 N/m jl7lll8l |. 

The resonant frequency of a mass-loaded cantilever is 
determined by the spring constant k and the mass of the 
test mass m t : 

cJo=fc/m t , (2) 

where too is the angular frequency of the first bending 
mode of the cantilever. The spring constant of each can- 
tilever was deduced from the resonant frequency, which 
was measured very precisely, and the mass of the test 
mass, discussed in the next section. Adding the test 
mass reduced the resonant frequency of the cantilever 
to ~ 300 Hz. As found from the resonant frequency, the 
addition of the test mass increased the spring constant 
to ~ 0.007 N/m; this increase was due to the shortened 
effective length of the cantilever once the test mass was 
attached. 

Cantilevers exhibit a Lorentzian transfer function be- 
tween driving force and amplitude; when driven at the 
resonant frequency fo by force F, the maximum displace- 
ment x at the center of mass of the test mass is 

x(fo) = F(f Q )Q/k, (3) 

where Q is the quality factor. The quality factors of 
cantilevers used in this experiment were found to be as 
high as 80000 in vacuum and at cryogenic temperatures. 

The energy from the thermal environment provides a 
constant series of random-phase impulses at all frequen- 
cies to the device. The resulting motion of the cantilever 
shows the Lorentzian transfer function, with a force spec- 
tral density Sf on resonance of 



S f = ^/4kk B T/Qiu , (4) 

where T is the cantilever temperature and kg is the 
Boltzmann constant. We used floppy (low spring con- 
stant), high quality factor cantilevers at low tempera- 
tures to reduce this limit. Force measurements were av- 
eraged as long as needed to see a signal above noise or as 
long as was practical. 



B. Test Mass 

To fabricate test masses, gold was deposited (using a 
thermal evaporator) into molds of plasma-etched silicon. 
After polishing of the top surface of evaporated gold, the 
silicon was dissolved to release the test masses. To make 
magnetic test masses, 1000 A of Ni was evaporated before 
the gold. The test masses were designed to be prisms 
(50 x 50 x 30) /im 3 in size; with these dimensions, a test 
mass could be affixed to a cantilever with one of its larger 
faces flush against the cantilever and aligned with the end 
of the cantilever. 



TABLE I: Dimensions of the test mass. 



Parameter 


Value 


Error 


Units 


Width 


51 


1.5 


/im 


Length 


51 


1.5 


/im 


Rectangular height 


32 


1.5 


fim 


Rounded volume 


7800 


3000 


/tm 3 


Missing volume from side face 


4000 


1500 


/tm 3 


Missing volume due to porosity 


2.5 


2.5 


% 


Density of gold 


19.3 





g/cm 3 


Total Mass of test mass 


1.64 


0.13 





A group of thirty test masses was weighed (in sets of 
five or six) with a microbalance to determine the aver- 
age mass. Because many of these test masses were im- 
perfect specimens, this measurement was only used as a 
guide. After weighing, test masses were examined under 
a scanning electron microscope (SEM) and etched by a 
focused ion beam (FIB) to determine more precisely the 
typical dimensions and porosity of the test masses. All 
the nonmagnetic test masses were fabricated together; 
those examined in the FIB were assumed to be repre- 
sentative of the entire batch, including the one used in 
Cooldown 04. The SEM showed the test masses to be 
slightly larger than intended (as confirmed by the mi- 
crobalance measurements), with a rounded face where 
the evaporated gold was polished. Etching with the FIB 
showed the evaporated gold to have some porosity near 
the side faces of the test mass. 

After data acquisition, the test mass and cantilever 
used in Cooldown 04 were examined under an SEM 
(Fig. • It was found that the test mass was mounted 
with its rounded side closer to the cantilever and with 
the top flat side tilted at 0.35 rad with respect to the 
cantilever in the y-z plane. Estimates of the size and 
uncertainty in the dimensions of the test mass were de- 
rived from examination of this test mass and of the larger 
group of masses from the same fabrication batch. SEM 
images of the Cooldown 04 test mass showed that the 
side face that was further away from the cantilever (due 
to the tilt of the test mass) was partially recessed; an 
estimate of this missing volume was included in the cal- 
culation of the mass. Error on these estimates was due 
to error in the SEM measurements and uncertainty in 
the exact shape of the rounded part on the polished side. 
The density of the gold was assumed to be the bulk den- 
sity of 19.3 g/cm 3 ; the observed porosity was included in 
volume and mass uncertainty. The dimensions and den- 
sity of the test mass and the experimental uncertainties 
in these parameters are shown in Tbl.[I] The mass given 
in Tbl. J] agrees well (within the given experimental un- 
certainties) with the mass deduced from a comparison of 
the resonant frequencies of the cantilever with the mass 
and neighboring cantilevers without masses. 



FIG. 3: Scanning electron micrograph of the test mass and 
cantilever used in Cooldown 04. A piece of dust is visible on 
one corner of the test mass. The cantilever is 50 /im wide. 



C. Drive Mass 

The drive mass was also fabricated by evaporat- 
ing gold into a mold of silicon. After polishing, the 
silicon substrate was diced into dies approximately 
1.8 mm x 1.3 mm in size. 

Gold and silicon have differing electrical conductivities, 
as well as differing mass densities. To eliminate the pos- 
sibility of a periodic coupling between the masses due to 
the Casimir force or charging of the drive mass, the pat- 
tern of the drive mass was buried beneath a thin ground 
plane. To bury the drive mass pattern, the polished side 
of the die was mounted on a quartz backplane, which 
became the bottom of the drive mass. The bulk of the 
remaining silicon on the top was removed, leaving a layer 
less than 2 /zm thick. On top of this layer of silicon was 
deposited a thin layer of aluminum oxide (for electrical 
insulation) and 1000 A of gold on top of an adhesion 
layer of titanium. This gold film was continued around 
to the side of the drive mass, where electrical contact was 
made to a ground wire glued on with silver epoxy. This 
thin ground plane masked variations in the Casimir force 
without notably affecting the varying gravitational field 
of the drive mass. 

As shown in Fig. 0] the main part of the drive mass 
pattern comprised five sets of gold and silicon bars, each 
1 mm long and 100 /im in each of the cross-sectional 
dimensions. The rest of the drive mass pattern pro- 
vided leads to which electrical contact could be made 
for grounding the meander (for the gravitational exper- 
iment) or for passing electrical current through the me- 
ander (for the magnetic experiment). 

A drive mass similar to the one used for measurement 
was etched by the FIB and examined under the SEM. 
On the polished side of the die, each gold bar had a band 
of indeterminate composition on either side, where the 
polishing created a wedge of silicon, gold, and polishing 
grit mixed together. This wedge extended 10-20 /jm into 
the gold bar at each gold/silicon boundary, tapering off 
at a depth of ~ 10 /im into the bar; this band of indeter- 
minate composition thus comprised a small part of the 
(100 x 100) /im 2 cross-section of each gold bar. After final 
fabrication and mounting of the buried drive mass, this 
polished side of the drive mass was facing away from the 



FIG. 4: [color online] Optical micrograph of the polished side 
of a drive mass showing the meander pattern. In the exper- 
iment, the polished side of the drive mass was facing away 
from the test mass. 

test mass. Because the Yukawa potentials in question are 
short-range, this imperfection at a distance of > 100 /im 
away from the test mass would have little effect on the 
results. In the analysis, this imperfection was taken as a 
small error on the bulk density of the gold in the drive 
mass. It was assumed that the drive mass had poros- 
ity on the sides of the gold bars similar to the amount 
observed in the etched test mass. 



D. Shield Membrane 

The test mass and cantilever were isolated from elec- 
trostatic and Casimir excitations by a shield membrane 
between the cantilever wafer die and the drive mass. The 
cantilever was held within a silicon wafer die approxi- 
mately 1 cm 2 in size. This die was glued to another 
silicon die, which was etched into a frame bearing a 3- 
/zm thick membrane of silicon nitride across an area of 
5.2 mm x 2.8 mm. The entire shield wafer die, including 
the membrane, was coated with gold on both sides, with 
a ground wire attached to one corner of the die. Due 
to t he g eometry and the tensile stress in the membrane 
[Tfll l2flL l2lj. the membrane was expected to be much 
stiffer than the cantilever. 

If there was some force between the drive mass and 
the shield at /o, the shield would move at Jq. The shield 
motion could drive the cantilever capacitively or via the 
Casimir force. However, because of the stiffness of the 
shield and the separation of the drive frequency from 
the signal frequency, any interaction between the shield 
and the drive mass would likely be too small to make the 
shield deflect enough to drive the cantilever a measurable 
amount on resonance. 



E. Piezoelectric Bimorph Actuator 

A piezoelectric bimorph actuator (hereafter referred to 
as the "bimorph" ) was used to move the drive mass Ion- 
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gitudinally underneath the test mass. The bimorph was 
driven at a subharmonic of the cantilever resonant fre- 
quency /o; the particular subharmonic (either /o/3 or 
/o/4) was chosen so the drive frequency was below the 
resonance of the bimorph but high enough to gain a res- 
onant enhancement in the amplitude and a reduction of 
nonlinearities in the bimorph motion. At low tempera- 
tures, self-heating of the bimorph, resonant enhancement 
of the motion, and driving voltages larger than the room 
temperature limits for the device allowed an amplitude 
of 100-125 jum of motion at a drive frequency f d of 90- 
120 Hz. Finite element analysis showed that the magni- 
tude of the time- varying force at 3/d from the Newtonian 
and any Yukawa potential between masses would vary as 
a function of the bimorph amplitude, with the maximum 
occurring at a bimorph amplitude of ~ 135 fim. 

The bimorph was secured at its base in a brass and 
Cirlex [2^ clamp. On one side of the bimorph was glued 
a capacitive electrode facing a counter-electrode mounted 
on the clamp. On the other side of the bimorph was 
mounted a small mirror. Prior to being mounted in the 
probe, the bimorph was calibrated. Room-temperature 
calibration of the bimorph was performed by measuring 
the rms capacitance between the two electrodes with a 
given driving voltage on the bimorph. The corresponding 
amount of motion was determined by reflecting a laser 
beam off the mirror onto a linear CCD array. At low 
temperatures, the amplitude of the bimorph motion was 
determined from a measurement of the rms capacitance. 

Uncertainties in the bimorph motion were tabulated 
through the linear fits used to compare the capacitance 
measurement in situ to the calibration. The total ex- 
perimental uncertainty in the amplitude of motion was 
10 /im, mostly a result of uncertainty in the original mea- 
surement of the laser spot position from the CCD array. 

A small ground cap was glued on top of the bimorph, 
to help isolate the shield membrane from the large ac 
voltages used to drive the bimorph. The drive mass was 
glued on top of this ground plane. 

For purposes of the calibration, the bimorph bending 
shape was considered to be an arc with constant radius 
along the length of the bimorph. For our analysis, it was 
assumed that the drive mass was moving purely in the 
horizontal plane; in fact as the bimorph bends, the drive 
mass will tip up at the ends, changing the vertical sep- 
aration between the test mass and the drive mass below 
it. Depending on the alignment of the masses, the re- 
sulting change in the vertical separation between masses 
may be as much as 2 /jm over the course of one period 
of drive mass motion. This effect was not expected to 
substantially change the signal and a full modeling of it 
has been left for future work. 



F. Vibration Isolation 

Measurement at the thermal noise limit required me- 
chanical excitation at the cantilever tip to be less than 



1 A (rms) on resonance; the mechanical excitation at the 
base of the cantilever was required to be reduced a factor 
of Q beyond this. In order to make a sensitive measure- 
ment, it was crucial to ensure that the bimorph was not 
simply mechanically exciting the cantilever. Even though 
the bimorph was moved at a subharmonic of the resonant 
frequency of the cantilever, nonlinearities in the bimorph 
motion (as in any piezoelectric device) meant that some 
component (typically a few percent) of its motion was 
at fa; this necessitated vibration isolation between the 
bimorph stage and the cantilever wafer stage. 

The bimorph and cantilever wafer were separated by 
two simple spring-mass vibration isolation stages. These 
stages each had resonant frequencies of ~ 2.5 Hz, cal- 
culated to yield roughly six orders of magnitude of at- 
tenuation at the bimorph frequency and eight orders of 
magnitude of attenuation at the signal frequency Jo- 

The electrical wires and the optical fiber that ran the 
length of the probe could have, if not loose enough, 
shorted out the vibration isolation system. Measure- 
ments of the cantilever at low temperatures proved to be 
the best test of the vibration isolation system. The me- 
chanical coupling between bimorph and cantilever could 
be assessed by measuring the motion of the cantilever 
with the bimorph moving at a large (~ 1 mm) vertical 
separation between masses, with the signal frequency on 
and off resonance of the cantilever. Such couplings were 
found to vary over the course of one cooldown, with no 
clear indication of failure of the vibration isolation. 

External vibration isolation protected the fragile parts 
of the probe from human mechanical disturbances when 
the drive mass was only microns away from the shield 
membrane. During data-acquisition, the entire cryogenic 
system was suspended from a thick concrete ceiling by 
springs, yielding a resonant frequency for the entire sys- 
tem of close to 2 Hz. 



G. Capacitive Position and Tilt Sensors 

Capacitive position and tilt sensors were used to indi- 
cate the relative position between the bimorph and wafer 
stages, so that the masses could be aligned to maximize 
the gravitational force between them. The capacitive po- 
sition sensor (CPS) was based on the design described 
in Ref. (23j. A quadrature pattern of gold electrodes 
was patterned on a quartz substrate and these electrodes 
were mounted on the bimorph stage. Above this, on the 
wafer stage, was mounted a single sense electrode. An 
ac voltage was applied to the quadrature electrodes and 
the current was read from the counter-electrode via a 
dual-phase lockin amplifier. With a phase shift of 90 
degrees between the signals applied to each of the four 
quadrature electrodes, the two channels of the lockin am- 
plifier provided readings linear in the relative x and y po- 
sitions between the bimorph stage and the wafer stage. 
With the signal applied in-phase to each of the quadra- 
ture electrodes, the signal provided an indication of l/z c , 
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where z c is the vertical separation between the sense and 
quadrature electrodes. 

The tilt sensors were simpler, each comprising an elec- 
trode mounted on the bimorph stage and a counter- 
electrode on the wafer stage above. In each case, an 
ac voltage was applied to one electrode, with the cur- 
rent read from the corresponding counter-electrode via 
a lockin amplifier. The tilt sensors provided an indica- 
tion of the vertical separation between the stages at two 
other locations in addition to the CPS. The geometry of 
the three capacitive sensors is shown in Fig. [S] 

Together, the capacitive sensors provided an indica- 
tion of the relative position and tilt between the stages, 
encompassing five degrees of freedom. The relative rota- 
tion between the stages about the z-axis was fixed during 
probe assembly. 

The capacitive sensors were calibrated at room temper- 
ature by raising the bimorph stage with a three-axis posi- 
tioner and recording the micrometer readings from the z- 
axis of the positioner along with the capacitive readings. 
The CPS was also calibrated in x and y using a similar 
technique. With the probe under vacuum, the three-axis 
positioning stage tilted due to atmospheric pressure; this 
tilt did change the correspondence between micrometer 
and CPS readings. Alignment between masses was es- 
tablished in room temperature at atmospheric pressure. 
At low temperatures, with the system under vacuum, the 
capacitive sensors and the room temperature calibrations 
(rather than the micrometer readings) were used to indi- 
cate the relative position and tilt between stages. 



H. Cryogenic Apparatus 

Cooling the cantilever was crucial to achieving a high 
force sensitivity. The probe was sealed in a vacuum can 
inside a 4 He cryostat, with an exchange gas space sep- 
arating the probe from the mechanical vibration of the 
boiling helium. The three-axis positioner was located 
outside of the cryostat, attaching to the bimorph stage 
via a vacuum feedthrough. The interferometer was also 
located outside the cryostat, with a length of optical fiber 
running down the probe. 

Base temperature of the probe was typically 10 K, 
though this temperature was raised a few degrees by the 
use of spring-adjustment heaters (described in Sec. lIVE|l . 
The noise temperature of the cantilever was found to be 
typically ~ 10 K more than the base temperature of the 
probe. 

Hold time for the cryostat was typically 4-6 days. Af- 
ter each helium transfer, all preliminary tests of the sys- 
tem were performed. Each data set presented in this 
paper was recorded entirely within one helium transfer 
of the respective cooldown. 




FIG. 5: Drawing of the bottom half of the probe. The bi- 
morph stage and wafer stage are shown with the wafer stage 
raised and rotated about the x-axis to reveal the geometry 
of the capacitive sensors. In practice, the two electrodes of 
a given sensor were separated by less than 1 mm in the z- 
direction. The three sensors were separated by ~ 4 cm in the 
x-y plane. The pads used for making electrical contact to the 
capacitive sensors are omitted from this drawing. 



I. Interferometer 

Cantilever motion was measured with a fiber interfer- 
ometer, based on the design described in Ref. |24|. An 
InGaAs laser diode sent hundreds of microwatts of 1310- 
nm light through a bidirectional fiber coupler, leaving 
1% of this light for the length of fiber that went to the 
cantilever. The cleaved end of the fiber in the probe was 
aligned to the test mass, forming a low-visibility Fabry- 
Perot cavity with a length of ~ 50 /im. Interference be- 
tween the light reflected off the test mass and the light 
reflected from the cleaved end of the fiber was measured 
with a photodiode via a transimpedance amplifier with a 
10 MO feedback resistor. 

Because of the low reflectivity of the cleaved end of 
the optical fiber, beam divergence, and imperfect align- 
ment between the fiber and the test mass, only one reflec- 
tion from each surface (the fiber end and the test mass) 
was expected to contribute to the interferometer signal. 
Thus, the dc interferometer signal was a sinusoidal func- 
tion of the distance between the fiber end and the test 
mass, modulo half the wavelength of the laser. The sen- 
sitivity of the interferometer was maximized when the 
distance between the cantilever and the optical fiber was 
adjusted to be at one of the points of maximum slope of 
this sinusoidal fringe (the center of the fringe). 

Cantilever motion was typically on the order of 
angstroms at low temperatures. For motion much less 
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than the wavelength of the laser, the interferometer sig- 
nal at a given frequency Vi (/) was linearly related to the 
amplitude of cantilever motion at that frequency x(f). 
The conversion between the voltage signal from the inter- 
ferometer and cantilever motion depended on the peak- 
to-peak amplitude V pp of the sinusoidal interferometer 
fringe and the relative position between the cantilever 
and the optical fiber that determined the location on this 
fringe of the interferometer dc level. When the distance 
between the fiber and the cantilever was adjusted to the 
center of the fringe, the amplitude of cantilever motion 
was determined by 



<f) = Vi{f) 



X, 



2irV, 



(5) 



where A; is the wavelength of the laser. 

Temperature control and high-frequency (> 100 MHz) 
modulation j2j| were used to stabilize the laser during 
data acquisition. The high-frequency modulation short- 
ened the coherence length of the laser, reducing the im- 
portance of stray reflections from connectors within the 
optical part of the interferometer circuit. 



J. Cantilever Position Adjustment and 
Characterization 

A piezoelectric stack (hereafter referred to as "piezo 
stack") mounted underneath the cantilever wafer was 
used to adjust the distance between the cantilever wafer 
and the optical fiber to maintain alignment at the center 
of the interferometer fringe. This piezo stack was also 
used to excite the cantilever for the purposes of charac- 
terization. 

In the probe, the piezo stack was mounted between the 
wafer frame and the wafer stage. Upon cooling, the stain- 
less steel wafer frame and stage would contract while the 
piezo stack would lengthen slightly. At low temperatures, 
it was expected that the differential thermal contraction 
of the materials would have made the wafer frame tilt in 
the x-z plane. A correction for this tilt was included in 
the room temperature alignment between masses. 



III. FINITE ELEMENT ANALYSIS AND 
SPATIAL PHASE-SENSITIVE DETECTION 

In order to deduce a bound on Yukawa-type deviations 
from Newtonian gravity, the measurement must be com- 
pared to the expected force between the masses for both 
a Newtonian potential and a Yukawa potential. Finite 
element analysis (FEA) with a mesh size of 5 /im was 
used to calculate the expected dc (without time varia- 
tion) gravitational force between the masses for a range 
of longitudinal (along the y-axis) positions between the 
masses, with the vertical separation held constant. 

At a given y-point (with other alignment parameters 
set), the dc Newtonian force Fj^ and dc Yukawa force 



i^ dc) (with a = 1) were calculated by sums over the 
drive mass and the test mass in the 5 /tin mesh. Only 
the vertical component of the force was considered and 
the attractive force was taken to have a positive sign. 
The Newtonian force sum is given by: 
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Here, the two sums are over the volumes of the drive mass 
and the test mass in the given mesh, Snid refers to the 
mass of the (5 /mi) 3 block in the drive mass, Suit refers 
to the mass of the (5 /urn) 3 block in the test mass, and 
r is the distance between the centers of mass of these 
two blocks in the summation. The vertical separation 
between the two mass blocks is z; the final term in the 
equation is a projection onto the vertical axis. Similarly, 
the Yukawa force (for a = 1) for a given value of A was 
calculated from the sum: 
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A time-variation (accounting for drive mass motion) 
was applied to this calculated set of dc forces as a func- 
tion of y-position and from this, the expected ac (finite 
frequency) force was extracted with a Fourier transform, 
yielding ac Newtonian and Yukawa forces. In the case of 
Cooldown 04, the third harmonic of the drive frequency 
was studied; the Fourier component of the measured force 
at this frequency is referred to as the third harmonic ac 
force. 

There were several inputs to the FEA model: the ge- 
ometry and density of the masses, the (x, y, z) position 
and tilts between masses, amplitude of bimorph motion, 
and the range A of the Yukawa potential being modeled. 
The output of the model was the ac Newtonian gravita- 
tional force between masses Fm and the ac Yukawa force 
Fy for a — 1 and a given A. For arbitrary a, the Yukawa 
force could simply be scaled by a. 

Both the magnitude and phase of this Fourier compo- 
nent with respect to the drive frequency are important. 
Clearly, the dc gravitational force between the masses 
reflects the periodicity of the drive mass pattern. In 
fact, the ac force also shows this periodicity. As shown 
in Fig. the magnitude of any ac gravitational force 
has a clear periodicity of 100 /im as a function of the 
//-equilibrium position of the drive mass with respect to 
the test mass, corresponding to the 100 /on half-period of 
the drive mass pattern. Each minimum of the third har- 
monic ac force magnitude is zero and is accompanied by 
a discontinuous phase change of n. The fourth harmonic 
ac force (not shown) has the same periodicity, though 
the fourth harmonic force has minima where the third 
harmonic force has maxima, and vice-versa. 

To exploit this geometric feature of our design, data 
were recorded at several values of the //-equilibrium po- 
sition between masses, scanning over more than 200 /im, 
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FIG. 6: The calculated magnitude (top) and phase (bottom) 
of the predicted third harmonic ac gravitational force between 
masses. The larger Yukawa force for a — 5 and A = 34 /im 
is shown by the dashed line; the smaller Newtonian force is 
the solid line. The phase is the same for a Newtonian or a 
Yukawa force. In this calculation, the drive mass was taken to 
be 27 /im underneath the test mass, without tilt. The middle 
gold bar of the drive mass is centered under the test mass at 
y-position of /nm. 

the period of the drive mass pattern. Comparison of 
these y-scan measurements to FEA predictions allowed 
us to discriminate between couplings that could be grav- 
itational in origin and spurious backgrounds that do not 
follow the gravitational pattern. By "locking-in" this way 
to the expected spatial periodicity, we were able to set a 
stronger and more accurate bound on a(X) than what a 
single force measurement would have provided. 



allowed selection of the better test masses from the fab- 
rication batch. 

After attaching test masses to cantilevers, the can- 
tilever wafer was glued into a stainless steel wafer frame. 
The shield wafer was then glued to this cantilever wafer, 
using a press to keep the wafers parallel. Measurements 
before and after each gluing showed the degree of tilt 
(<C 1 mrad) between the two wafers. Photographs of the 
two wafers before and after gluing were used to deter- 
mine the relative position between each cantilever on the 
wafer and the center of the shield membrane. 

By design, the shield membrane was 10 /im below the 
surface of the shield wafer that was glued to the cantilever 
wafer. The glue between the wafers added typically 1- 
5 /im to this distance. The offset of the shield and the 
thickness of the glue layer did limit the minimum vertical 
separation between masses. However, reducing this dis- 
tance significantly would cause the cantilever to snap-in 
and adhere to the shield [2(| . 

B. Probe Assembly 

Before mounting in the probe, the shield membrane 
was carefully examined for dust. Any pieces of dust with 
a resolvable height profile (> 2 /im) were removed with 
a microprobe. Then the wafer was mounted in the wafer 
stage and the optical fiber was aligned to the test mass. 
Calibrations of position sensors were performed and the 
drive mass and the bimorph were mounted in the probe. 

C. Drive Mass Gluing 



IV. EXPERIMENTAL METHODS 

In order to accurately compare measurements to the 
FEA, alignment between masses had to be known as ac- 
curately and precisely as possible. In order to maximize 
the gravitational force between masses, the goal was to 
center the masses with respect to each other with no tilts 
about the y-axis (0 XZ ) or the x-axis {0 yz ). The align- 
ment coordinates are diagrammed in Fig. {7\ With no 
tilts in the system, the third harmonic ac gravitational 
force would be maximized with the test mass at x = /im 
and y = 50 /im with respect to the center of the drive 
mass at equilibrium. 

Alignment between the masses was fixed by room- 
temperature preparations and low-temperature adjust- 
ments. This process began with the assembly of the can- 
tilever and shield wafers. 



A. Wafer Assembly 

The test mass was attached to the cantilever using 
epoxy applied with a microprobe under an optical mi- 
croscope. The optical microscope used for this process 



The drive mass was glued to the bimorph with the 
bimorph positioned so that the drive mass was pushed 
against the silicon shoulder of a shield wafer, mounted in 
the probe. This process, only repeated when the bimorph 
or the drive mass needed to be replaced, achieved approx- 
imate parallelism between the drive mass and any can- 
tilever wafer. Before each cooldown, parallelism between 
the shield wafer and the drive mass was again examined 
and adjusted. Uncertainties in the tilt between the drive 
mass and the cantilever were dominated by uncertainties 
in this process. 

D. Alignment 

After the drive mass was glued to the bimorph, the 
parallelism between the shield wafer and the drive mass 
was examined optically, using images recorded from a 
telescope via a CCD camera. Parallelism was assessed in 
the x-z and y-z planes by looking at the reflection of the 
drive mass in the gold-coated shoulder of the shield wafer. 
CPS readings and optical images were compared as the 
bimorph stage was moved to change the vertical separa- 
tion between the drive mass and the shield wafer. Tilts 
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FIG. 7: [color online] Schematic showing the geometry and 
alignment parameters between the test and drive masses. On 
the left, the top view shows the x-y plane; the origin, marked 
with an asterisk, is the center of the drive mass in this plane. 
With the drive mass at equilibrium, the x-y position of the 
test mass was defined with respect to this origin. This fig- 
ure shows the masses to scale; the shield membrane and 
the ground plane over the drive mass are omitted. On the 
right, the side view shows (not to scale) the z-x plane. The 
z-separation between masses is the face-to-face vertical dis- 
tance. The test mass tilt, exaggerated in this figure for the 
purposes of illustration, is the tilt of the test mass with re- 
spect to the cantilever. The angle between the drive mass and 
the cantilever in this plane was defined as 9 XZ . 

measured in two vertical planes (9 XZ and 9 yz ) were ad- 
justed using turnbuckles on the lower vibration isolation 
stage (the wafer stage), with additional compensation for 
the expected piezo-induced tilt of the wafer frame upon 
cooling. Uncertainty in this part of the alignment was 
dominated by optical limitations. 

The large three-axis positioner was used to adjust the 
alignment of the drive mass with respect to the test 
mass both at room temperature and at low temperatures. 
Micrometers on this stage were accurate to 2.54 fxm 
(0.0001 in), with the y and z micrometers controlled by 
motors. 

After optimization of the tilt, the bimorph frame was 
moved so that the drive mass was centered in the x-y 
plane with respect to the shield membrane; a telescope 
with a reticle aided this centering process. The bimorph 
stage was then slowly raised until the drive mass con- 
tacted the shield. This contact gave an impulse to the 
shield and the cantilever, creating a clear signal on the in- 
terferometer. CPS and tilt sensor readings were recorded 
several times during this process; these alignment points 
were the targets for initial alignment when the system 
was cold. Uncertainties in this process included optical 
limitations and motion of the three-axis positioning stage 
in the x-y plane as it was raised along the z-axis. 

The drive mass and shield were examined through a 
telescope for any indication of dust on either surface be- 
fore the vacuum can was closed. 



E. Re-Alignment 

Under vacuum, the tilt of the positioning stage con- 
nected to the bimorph frame significantly changed the 



alignment between the wafer stage and bimorph stage, 
and hence the alignment between masses. At base tem- 
perature, the tilt of the wafer stage was adjusted by heat- 
ing two of the three springs on the upper vibration iso- 
lation stage. To heat a spring, a current on the order of 
20 mA was passed through a manganin wire (of resistance 
~ 50 SI) wrapped around the coils of the BeCu spring; 
such heating could increase the length of the spring on 
the order of 100 fim. Temperature and position of the 
probe stabilized after several hours. 

The long time scale of the spring heating prohibited 
fine adjustment of the tilt. However, the large separa- 
tion between the tilt sensors (~ 38 mm) in the x-y plane 
in comparison to the size of the drive mass (~ 1.5 mm) 
allowed for a coarse readjustment of tilt at low temper- 
atures; typically, the vertical separation between the tilt 
sensors for a given CPS z c -reading was within 100 /Ltm of 
the room temperature alignment points. 

With the tilt adjusted, the x-y position of the bimorph 
stage was adjusted to regain the room temperature align- 
ment points. The manual operation of the ^-micrometer 
necessitated coarser alignment in this direction, since the 
micrometer could not safely be adjusted with the drive 
mass positioned close to the shield wafer. However, the 
alignment between the masses in the ^-direction only 
needed to be accurate to within a couple of hundred mi- 
crons, since the Yukawa forces being studied were short 
range in comparison to the 1-mm length of the drive mass 
bars. The motor control of the y and z micrometers al- 
lowed adjustments in these directions to within a micron, 
with uncertainty being dominated by noise (electrical and 
mechanical) on the capacitive readings. 

Additional information about the alignment between 
masses was provided by moving the bimorph frame in 
the y-direction until the drive mass contacted the side 
of the silicon frame bearing the shield membrane. This 
test was performed with the bimorph static and lowered 
so that the drive mass was ~ 200 fim away from the 
shield. The contact between the drive mass and the wafer 
gave a mechanical impulse to the cantilever, clearly vis- 
ible on the interferometer signal. The distance between 
the alignment point and this contact point provided ad- 
ditional confirmation of the drive mass alignment along 
the y-axis. 

After realignment, the bimorph frame could then be 
moved to the y-position, determined by the photographs 
of the wafer, at which the drive mass was centered un- 
derneath the test mass. The experimental values and 
uncertainties in the alignment parameters are given in 
Tbl.|H| 



F. Vertical Positioning of Drive Mass 

Typically, gravitational measurements were recorded 
with the drive mass positioned 10-15 /jm away from the 
shield, with this safety factor allowing for the small ver- 
tical motion of the drive mass over the course of the bi- 
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TABLE II: Alignment between drive mass and test mass 



Parameter 



Value Error Units 



s-position 
y-position 
z-separation 
\0 X x \ x width 6 
\8 yz \ x length' 



-91 

(-12)-(283) 
24-28 
1 
6 



110 

119" 
3.4 
6 
9 



/im 

fim 
fim 



"This is the error on the location of the entire range. Within this 
range, there was an uncertainty of about 5 fim on each point. 
''Here, width is of the drive mass die: 1.3 mm. 
c Here, length is of the drive mass die: 1.8 mm. 



morph swing, backlash in the motor used to operate the 
z-axis of the three-axis positioner, and the possibility of 
drift during a data run. 

To determine the vertical (z-) separation between the 
masses, the bimorph frame was slowly raised until the 
drive mass contacted the shield, giving an impulse to 
the cantilever clearly visible on the interferometer. The 
bimorph was then lowered to the desired distance from 
the shield with the CPS indicating the amount of motion. 
The largest sources of uncertainty in the z-separation 
between masses were bouncing during contact and data 
acquisition and the possibility of dust on the shield or 
drive mass. Due to lack of precision in setting a given 
z-separation (due to limitations of the motor driving the 
z-micrometer) , the z-separation during data acquisition 
varied over the course of a day; this is the reason for the 
range given in Tbl. [H] 



G. Cantilever Characterization 

After reaching base temperature, the piezo stack was 
used to verify the interferometer fringe. The shape of the 
fringe was an important indicator of alignment between 
the optical fiber and the test mass. In earlier designs 
of this apparatus, the fiber alignment often drifted upon 
cooling so that more light was reflecting from the shield 
than from the test mass. Other times, light reflected 
from the cantilever as well as from the test mass. The 
former situation created a large background level on the 
interferometer; the latter would make it impossible to 
record data. In Cooldown 04, the shape of the fringe 
indicated no misalignment between the optical fiber and 
the test mass at low temperatures. 

The resonant frequency fo was determined by exciting 
the cantilever at large amplitudes with the piezo stack 
and comparing the amplitude and phase of the resulting 
signal to the drive signal. The quality factor was de- 
termined by timing the ringdown of the cantilever from 
excitation on resonance. The fringe height, the resonant 
frequency, the quality factor, and the dc level of the in- 
terferometer on the fringe were checked between data col- 
lection runs. 



H. Bimorph Actuation 

After alignment between the masses was established, 
the bimorph was lowered to a safe distance and an ac 
voltage was applied to the bimorph via a high voltage 
amplifier. The capacitive measurement of the bimorph 
motion indicated that the bimorph motion for a given 
driving voltage increased over time, typically coming to 
equilibrium within 30 min. 



I. Experimental Degrees of Freedom 

To make the most sensitive gravitational measurement, 
data were recorded as a function of y-position between 
masses, with a small vertical separation between masses 
and the drive frequency fd tuned to be /o/3. Diagnostics 
included data runs with fd tuned slightly off this reso- 
nance, with the drive mass far away from the test mass, 
with the bimorph not moving, or with fd = jo/4. 



V. DATA ACQUISITION AND AVERAGING 

Two streams of time-series data were recorded via an 
analog-to-digital converter (ADC) on a personal com- 
puter: the voltage from the interferometer and the 
voltage from the function generator that was driving 
(through a high-voltage amplifier) the bimorph. The in- 
terferometer signal showed the cantilever motion. The 
function generator signal (at frequency fd close to 
100 Hz) provided an important timing signal. Data were 
recorded at a frequency of 10 kHz. Before the ADC, the 
interferometer signal was ac-coupled to a pre-amplifier 
with a high-frequency rolloff at 3 kHz to avoid aliasing. 

Any motion of the cantilever that was due to the mov- 
ing drive mass would be at a definite phase with respect 
to the drive mass motion and at a harmonic of the drive 
frequency. The function generator signal (the drive sig- 
nal) provided a proxy for the drive mass motion and the 
analysis of the cantilever motion then included a phase 
defined with respect to this drive signal. 

Each data run included data recorded over a period of 
time t tl with time-series data examined in shorter records 
of time to. The time to, typically on the order of Q/ fo, 
was chosen so that each record of length to could be con- 
sidered statistically independent of other records, while 
maximizing the total number of records N = it/to- 

For each record of length to > the time-series data were 
truncated to include an integer number of periods of the 
function generator signal. This truncation shortened the 
length of each record a negligible amount and increased 
the accuracy of the amplitude and phase of the harmonics 
of fd reported by a Fourier-transform of the record. 

The drive signal, resulting from a very clean function 
generator signal, showed an unambiguous peak on the 
Fourier transform at the drive frequency fd- The har- 
monic of interest of this drive frequency was then selected 
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from the Fourier-transform of the interferometer signal. 
In the case of Cooldown 04, measurements were made 
from the third harmonic of the drive signal, intended to 
be equal to /o. The signal from the interferometer was 
converted to an amplitude of cantilever motion, as de- 
scribed in Eq. From this, the third harmonic ac force 
on the cantilever (assuming the 3/d = /o) was determined 
using Eq. [3J 

The magnitude and phase (or, equivalently, the real 
and imaginary parts) of the Fourier transform are im- 
portant to consider. For each record t , the real Fr and 
imaginary Fj parts of the Fourier transform at 3/^ were 
found, providing N measurements of Fr and Fj. Aver- 
aging each of these two sets provided the means Fr and 
-Fj for the given data run, with a standard error on each 
mean determined by the standard deviation of the set, 
reduced by a factor of y/N — 1. 

Aside from the statistical error, uncertainty in the as- 
sessment of the force on the cantilever from this averaging 
technique resulted from uncertainty in the fringe height 
of the interferometer, the interferometer dc position rela- 
tive to the center of the fringe, the cantilever Q, and the 
setting of fa to be exactly /o/3. Furthermore, Eq. as- 
sumes that the measured displacement of the cantilever 
is at the center of mass of the test mass. Due to the 
mode shape of the cantilever, described analytically in 
Ref. j2]} and modeled in ANSYS, the deflection at either 
end of the test mass would be different by ~ 15-20% from 
the deflection at the center of mass. Though the qual- 
ity of the interferometer signal indicated that the optical 
fiber was focused on the test mass (instead of the can- 
tilever), the degree of freedom in the fiber position along 
the length of the test mass provided another uncertainty 
in the conversion between interferometer signal and force. 
These uncertainties are summarized in Tbl. IIIII 



A. Feedback Cooling 

A high quality factor was essential to achieving a sen- 
sitive force measurement. However, a high quality factor 
also made measurement difficult; the resonant frequency 
had to be determined to an accuracy of ~ fa/('SQ) and 
between data runs or after any excitation of the can- 
tilever, a waiting period of at least ~ 3Q//o had to be ob- 
served to allow the cantilever to ringdown. In Cooldown 
04, Q ~ 80000. To facilitate the data-acquisition process 
while maintaining the low thermal noise of the high qual- 
ity factor cantilever, feedback cooling was used to reduce 
the effective quality factor and the noise temperature T 
of the cantilever in its lowest mode. 

The feedback-cooling apparatus was based on the de- 
signs of the Rugar group |28| and the feedback circuit de- 
scribed in Ref. [jjj] . The interferometer signal was phase- 
shifted, attenuated, and used to drive the piezo stack. In 
this way, the cantilever was driven with a phase-shifted 
version of its own thermal noise. The circuit was adjusted 
to incorporate a phase shift of n over the thermal noise 



TABLE III: Uncertainties in measured force magnitude 



Parameter 


Value" 


Error Units dF b 


(%) 


Fringe height 


0.228 


0.002 


V 


1 


Fringe position from center 





10 


% 


2 


/o c 


331.178 


0.002 


Hz 


2 


Q 


10000 


1000 




10 


A: 


0.0071 


0.001 


N/m 


8 


Fiber alignment 





25 


)im 


10 


Total 






%N 


17 



"Values given are typical for Cooldown 04, with feedback cooling 
on. For parameters (such a fringe height) that varied over the 
course of a day, the error given does not include long time scale 
drift, only the uncertainty over the course of an hour-long data 
run. 

''This is the proportional error in the force measurement resulting 
from the error in the named parameter. 

c The amplitude of the Lorentzian transfer function is extremely 
nonlinear as a function of frequency. Because the error considered 
here makes an assumption of linearity, the error in /o is overesti- 
mated to account for the nonlinearity; the relative error in the force 
was taken to be half of the change in the measured force that would 
result from twice the experimental uncertainty in the frequency. 



bandwidth of the first mode of the cantilever. Over this 
small bandwidth, this additional phase shift turned the 
feedback into negative feedback on the velocity (rather 
than on the position) of the cantilever. A cantilever is a 
damped simple harmonic oscillator; negative feedback on 
its velocity increases the damping, decreasing Q. It also 
decreases the amplitude of all thermal noise excitation 
in this first mode in a way that decreases the effective 
temperature of the cantilever in that mode. The thermal 
noise spectra of the cantilever in Cooldown 04 with and 
without feedback are shown in Fig. |SJ 

Within the accuracy of the measurement of these two 
parameters, the ratio of Q/T was maintained with this 
feedback and thus the low thermal noise limit was not 
affected by the reduction of the quality factor. However, 
the sensitivity of the cantilever, as described by Eq. |3 
does depend on Q; feedback cooling in Cooldown 04 did 
reduce the voltage signal-to-noise ratio of the measure- 
ment. 



VI. TESTS OF THE APPARATUS 

A. Thermal Noise 

To verify the functioning of the cantilever, to confirm 
the conversion from interferometer signal to amplitude 
of cantilever motion, and to determine the force limit 
in the experiment, thermal noise of the cantilever was 
recorded. During thermal noise measurements, the bi- 
morph was turned off. Thermal noise may be examined 
several ways. The noise temperature of the cantilever 
was determined from power spectra averaged over several 
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(F R + F/-) 1 / 2 due to thermal noise will decrease with 
the square root of the averaging time, as seen in Eq. 0] 

Even in the case of a coherent driving force on the 
cantilever, where the measured force at some phase is 
statistically distinguishable from zero, thermal noise is 
important to consider — thermal noise provides the statis- 
tical uncertainty on any force measurement. Averaging 
longer reduces the statistical uncertainty due to thermal 
noise. 



FIG. 8: [color online] Measured thermal noise data without 
(a) and with (b) feedback are shown in comparison to the 
Lorentzian curves derived from the measured k, Q, effec- 
tive T, and /n. The rms amplitude of cantilever motion is 
plotted as a function of frequency. The offsets added to the 
Lorentzian curves are judged from the power spectrum. The 
background in the figure on the right is higher because of the 
larger bandwidth. With feedback, Q ~ 10400 and the spec- 
trum is the result of averaging thirty power spectra from 2 min 
data records; the effective temperature was ~ 3.5 K. Without 
feedback, Q ~ 75000 and the spectrum is an average of thirty 
records of length 4 min each; the effective temperature was 
~ 30 K. 



consecutive time records. The equipartition theorem re- 
lates the mean squared displacement (a: 2 ) in one mode of 
the cantilever with spring constant k to the temperature 
T: 



(x 2 ) = k B T/k. 



(8) 



By summing the squared amplitudes of cantilever mo- 
tion over the thermal noise peak, the effective temper- 
ature was found. The white noise background of the 
interferometer was subtracted from this summation and 
the sum only included the thermal noise bandwidth so 
as to exclude electronic noise on the interferometer. The 
power spectrum of the cantilever due to thermal noise 
compared well to the expected Lorentzian, as shown in 
Fig- El Though the temperature of the cantilever without 
feedback was found to be ~ 30 K, higher than the tem- 
perature of the probe, the agreement between the mea- 
sured noise spectrum and the Lorentzian function using 
the measured values of Q, fc, and T confirmed our assess- 
ment of these important parameters. 

Thermal noise is random phase noise. Though a can- 
tilever is constantly excited by the energy in its thermal 
environment, this excitation comes as a series of random- 
phase kicks, with the ringdown from each kick character- 
ized by the quality factor. When the real and imaginary 
parts of a given Fourier component of thermal noise are 
considered, the measurements (over many time records) 
show a force that is statistically indistinguishable from 
zero. In this case, the mean of the set of measurements 
of Fr should be less than twice the standard error on 
the mean: Fr < 2an/y / N — 1, where an is the stan- 
dard deviation in the set of N measurements. The same 
is true for the set of measurements of Fj or any other 
phase of any Fourier component. The measured force 



B. Magnetic Analog Experiment 

In Cooldown 02, the buried drive mass and a mag- 
netic test mass were used in order to test the couplings 
between masses via a large, measurable magnetic force. 
The probe design was an earlier version than the one 
described above (used in Cooldown 04) and alignment 
methods were less precise. Nonetheless, this experiment 
served as an important demonstration of the functioning 
of the apparatus. 

To couple the drive mass to the magnetic test mass, 
an ac voltage could be applied to the drive mass with 
the bimorph not moving or a dc current could be drawn 
through the drive mass with the bimorph moving. In 
either case, the magnetic coupling between the drive 
mass and the test mass would reflect the 200 /j,m half- 
periodicity of the magnetic field across the drive mass. 
The exact direction of the dipole moment of the test 
mass, not known, would strongly influence both the mag- 
nitude of the coupling as well as the exact dependence 
of the force magnitude on the y-position between the 
masses. However, FEA showed that any magnetic force 
would have a distinct periodicity of 200 /im in the magni- 
tude of the force, with each force minimum accompanied 
by a phase change of n. These features are qualitatively 
similar to those of the gravitational force shown in Fig.El 
though the spatial period across the drive mass of this 
magnetic coupling would be twice that of a gravitational 
coupling. 

To quickly confirm the position of the drive mass with 
respect to the test mass, 0.5 V (~ 2 mA) at /o was ap- 
plied across the drive mass meander, with the bimorph 
not moving. With the drive mass close to the test mass, 
this created a large excitation of the cantilever, clearly 
visible on the interferometer signal viewed on an oscillo- 
scope. The drive mass was moved underneath the test 
mass and from this, the locations (in the CPS coordi- 
nates) of the distinct phase changes in the magnetic force 
were noted and used as reference in the measurements 
that followed. 

A //-scan was performed with the bimorph moving 
~ 100 //m in amplitude at /o/3. Measurements were 
made at intervals of ~ 25 /im in y, covering the entire 
drive mass pattern. At each y-position, measurements 
were made both with a current (0.5 mA) across the drive 
mass and with no drive mass current. Even with short 
averaging times (20 s with a current on and 10 min with 
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no current), a clear force was measured in each case. The 
results of this y-scan are shown in Fig. [5] 

With the current drawn across the drive mass, the force 
magnitude showed five distinct periods, corresponding to 
the five sets of gold and silicon bars. As expected, the 
phase of the force as a function of y-position showed dis- 
tinct discontinuous phase changes of it associated with 
each force minimum. Without knowing the direction and 
size of the magnetic dipole moment of the test mass, 
this measurement could not be used as a precise cali- 
bration. Furthermore, a broken shield prevented a pre- 
cise assessment of the vertical separation between masses. 
Nonetheless, this measurement provided a confirmation 
of the alignment between masses, the FEA predictions, 
and the relative position between data points as deter- 
mined by the CPS. 

Measurements with no current across the drive mass 
also showed a force with a distinct periodicity, though 
in this case the periodicity was 100 /im in the magni- 
tude of the force, as expected for a gravitational force. 
Phase changes of tt were also evident in this measure- 
ment, though the phase changes were continuous. Con- 
tinuous phase changes would be expected if there was a 
constant force in addition to the varying force across the 
y-scan. 

The measured force in the case of no drive mass current 
was much larger than any expected gravitational force . 
The shield and the buried drive mass would not permit 
any electrostatic or Casimir force to drive the cantilever 
in a way that would show this periodicity of the drive 
mass. This force could only be magnetic in origin. 

Silicon and gold have small diamagnetic susceptibili- 
ties, on the order of —0.16 x 10~ 8 m 3 /kg. Due to the 
differing mass densities of the gold and the silicon, each 
bar of silicon has a susceptibility that is smaller than the 
gold bar of the same size. Magnetized in any ambient 
field, the drive mass will have a spatially-varying mag- 
netic field even in the absence of applied current. This 
varying magnetic field will couple to the test mass and 
in this way, the magnetic test mass will act as a sus- 
ceptometer. Though this effect is small, our experiment 
is sensitive enough to measure such a coupling; a rough 
calculation shows that even in the earth's magnetic field, 
this effect can be as large as the measured force shown 
in Fig.HJ 

The functioning of the magnetic test mass as a sus- 
ceptometer limited the sensitivity of any gravitational 
force measurement with the magnetic test mass. How- 
ever, the measurement of both forces with the magnetic 
test mass — including both the 200 [im and the 100 ^tm 
periodicities in the force magnitude — provided an excel- 
lent verification of the sensitivity of the cantilever as a 
force detector and demonstrated that the experimental 
apparatus can measure the coupling between masses. 
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FIG. 9: [color online] The third harmonic force magnitude 
(a) and phase (b) as a function of y-position (relative to the 
first point) between the magnetic test mass and the drive 
mass with current on and off. The magnitude here is calcu- 
lated from the mean of the real and imaginary components: 

/ 2 2 
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C. Possible Background Couplings 

With a non-magnetic test mass, the buried drive mass 
and the shield eliminated the possibility of magnetic, 
electrostatic, or Casimir couplings to the test mass that 
could be mistaken for gravity. Couplings due to induced 
currents in both masses, magnetic impurities in the drive 
mass, or induced moments in both masses were all es- 
timated to be much less than the target thermal noise 
limit of 10 -18 N. The possibility of a coherent gravita- 
tional excitation of the cantilever due to anything other 
than the drive mass was made unlikely by the extremely 
high quality factor of the cantilever and the relatively 
high resonant frequency as compared to the motion of 
commonplace objects or people in the laboratory. 

If the vibration isolation was shorted, the bimorph 
could have shaken the cantilever. Measuring with the 
third harmonic 3/^ on- and off-resonance gave an indica- 
tion of whether the vibration isolation was failing. Can- 
tilever motion due to a mechanical excitation would be 
reduced by a factor of Q off-resonance. 

Finally, there could have been electrical or mechani- 
cal coupling between the bimorph and the interferome- 
ter that was unrelated to the cantilever. Such electri- 
cal or mechanical couplings to the interferometer would 
not vary as a function of the drive mass equilibrium y- 
position and thus would never be mistakenly interpreted 
as a gravitational signal. However, the resulting voltage 
noise could impede the detection of a small gravitational 
signal. 

Electrical coupling could result from imperfect ground- 
ing of the function generator or the circuitry of the laser 
or the interferometer. Mechanical coupling could have 
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resulted from the optical fiber being shaken by the mov- 
ing bimorph. In both cases, the coupling would produce 
a voltage noise on the interferometer that would be in- 
sensitive to small changes in frequency; the signal on the 
interferometer (with the subtraction of thermal noise) 
would be the same on- and off-resonance. Examinations 
of signals at fd and 2fd also provided a helpful diagnostic 
of such couplings. 

The optical fiber ran down the length of the tube to 
which the bimorph frame was attached. Motion of the 
bimorph could couple to the optical fiber itself, shaking 
the fiber at fd- Nonlinearities in the bimorph itself and 
in any part of the mechanical path between the bimorph 
and the optical fiber could result in the fiber also being 
shaken at harmonics of fd- Because there were reflections 
within the optical part of the interferometer circuit (due 
to the finite reflectivity of all fiber connectors) that cre- 
ated stray interferometric paths, vibration of the optical 
fiber could create a measurable signal on the interferom- 
eter at the frequency of vibration. Modulating the laser 
at high frequencies to reduce the coherence length could 
reduce but not eliminate this kind of noise. Having a 
large fringe height and a high quality factor to increase 
the sensitivity of the interferometer as described by Eq.[3| 
and Eq. [3] would reduce the importance of this kind of 
background noise. 



VII. EXPERIMENTAL RESULTS 

In Cooldown 04, a sensitive gravitational measurement 
was made using a nonmagnetic test mass, a buried drive 
mass, a new probe design to improve fiber alignment to 
the test mass, and more precise alignment techniques be- 
tween the masses. A y-scan covering almost 300 /im was 
recorded over two days. On each day, measurements were 
made at six points. The first data run on Day 1 was 48 
min long; all other data runs were 60 min in length. Data 
were recorded in records of 4 min in length, though analy- 
sis studied the records in 30 s segments. Feedback cooling 
was used to reduce the Q of the cantilever from ~ 80000 
to ~ 10000. The bimorph was oscillated at an amplitude 
(at the top surface of the drive mass) of 125 /im. 

The measured force and phase for the y-scan are shown 
in Fig. ^1 along with the thermal noise limit. On both 
days, the measured magnitude of the force showed an 
apparent periodicity of 100 /im, though there was no clear 
periodicity in the measured phase. All measured forces 
were close to thermal noise in magnitude. 

On Day 2, the measured force was consistently higher 
than on Day 1. The interferometer signals at fd and 2fd 
were also larger on Day 2, suggesting the presence of some 
mechanical or electrical coupling to the interferometer. 

The error bars in Fig. El show twice the statistical 
error on the measurements; the local maxima in the force 
measurements are statistically distinguishable from zero. 
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FIG. 10: [color online] Force magnitude (top) and phase (bot- 
tom) vs. y-position from Cooldown 04. Day 1 data are on the 
left and Day 2 data are on the right. The black dashed line 
shows the approximate level of measured thermal noise, which 
was slightly less than the predicted thermal noise level. Plot- 
ted points are the means for the full averaging time (60 min, 
except for the first measurement, which was 48 min); the error 
bars show twice the standard (statistical) error on the mean. 
Lines connecting the measured points are guides to the eye. 

VIII. ANALYSIS 

The measurements from Cooldown 04 do not clearly 
look like a gravitational force. Even if part of the mea- 
sured force was gravitational in origin, some of the mea- 
sured force was certainly due to thermal noise and pos- 
sibly due to some other background. To set a bound on 
Yukawa-type deviations from Newtonian gravity, these 
measurements as a function of y-position were compared 
to FEA predictions of the gravitational couplings be- 
tween the masses with the experimental conditions of 
Cooldown 04. 

A least-squares fitting, with a as a free scaling param- 
eter, of the predicted force as a function of y-position 
to the measurements provided a best-fit a for a given A. 
The best-fit a is based on the model that the measured 
force F m (y) is related to the Newtonian force F N (y) and 
the Yukawa force Fy{y) by 

F m (y) = F N (y) + aF Y (y, A) + F , (9) 

where y is the y-position between masses and F is some 
constant background. This equation was applied to the 
real and imaginary components of the force. MATLAB's 
"fminsearch" minimization routine was used for the least- 
squares fitting, with initial conditions chosen from a pre- 
liminary coarse-grained search for the fit parameters that 
would minimize the square error for a typical FEA result. 

This least squares analysis was performed separately 
on each day's worth of data. The real and imaginary com- 
ponents (rather than the magnitude and phase) of the 
force were considered in the fitting; the algorithm min- 
imized the summed square difference between the mea- 
sured force (considering both Fn(y) and Fj(y)) and the 
predicted force. Though each data run had 120 measure- 
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FIG. 11: [color online] Shown are the real and imaginary com- 
ponents of the Fourier transform at the frequency of inter- 
est for the second point of the Day 2 data — this is at the 
y-position at which the maximum force was measured for the 
day. Data here are reported in terms of the rms force, Fr and 
Fr. The small blue crosses each show the result from one 
30 s segment of time series data. Vertical error bars on two of 
these points indicate the statistical error, which is taken to be 
the same on each point and in all directions in this complex 
plane. The circle shows the mean along the two axes; the 
error bars on the mean are smaller than the size of the circle. 



merits each of Fr and Fi (except the first run, which 
had 96), the fit was performed using the means Fr and 
Fi at each y-point, in order to reduce computing time. 
Tests with real and simulated data showed that the re- 
sults were the same when the fit considered the entire set 
of data as when only the means were used. The scatter 
and the means of Fr and Fi at one y-point are shown in 
Fig. EU 

In addition to a, four other parameters were included 
in the fit. An offset y Q in the y-position accounted for 
the large uncertainty in the location of the range of mea- 
surements. Though the CPS gave an accurate indica- 
tion of the y-distance between two measurement points, 
there was an uncertainty of ~ 100 /xm (Tbl. [H} in where 
along the meander pattern the drive mass was located 
underneath the test mass. An offset in the phase 9q 
was included to account for the unknown relationship 
between the drive signal (from the function generator) 
and the phase of the signal measured on the interferom- 
eter. Though FEA predictions show how the phase of a 
measured coupling would vary, phase offsets due to cir- 
cuitry were not measured and thus the parameter 9o was 
necessary in the fit. 

Finally, an offset i?o on the real part of the force and an 
offset Iq on the imaginary part of the force were included 
in the fit. These two parameters account for the pos- 
sible presence of a constant force measured in addition 
to any varying force (which could account for the lack 
of a clear periodicity in the phase, even if a measurable 
gravitational force was present). Tests with simulated 
data showed that the inclusion of these additional offset 
parameters did not change the other best-fit parameters 



by unacceptable amounts. Moreover, on the average, the 
best fit i?o and Iq would be statistically indistinguish- 
able from if no offset was included in the simulated 
data. The Akaike Information Criterion [2^ indicated 
that using the additional offset parameters in the fit did 
in fact make a better model. Moreover, visual compari- 
son of the best fit FEA curves to the data indicated that 
inclusion of the offset improved the fit. 



A. Experimental Uncertainty 

There were many experimental uncertainties to be 
counted in determining the best-fit a(A). Many parame- 
ters, such as the geometry and density of the test masses, 
the amount of bimorph motion, and the vertical separa- 
tion and tilt between the masses, entered into the deter- 
mination of a only through the FEA and may have highly 
nonlinear effects on the best-fit a. Ideally, all uncertain- 
ties would be considered in a Monte-Carlo simulation. 
To make the computing tractable, only z-separation and 
tilts 9 XZ and 9 yz were varied in a Monte Carlo fashion. 
Other uncertainties were included after the least-squares 
results, as discussed below. Because our primary goal 
was to set an upper bound on deviations from Newtonian 
gravity, when simplifications were required, we chose to 
use "worst-case" estimates if possible, including errors in 
a way that would make the best-fit a(X) larger. 



B. Monte Carlo 

The FEA simulation was run 320 times, with the input 
parameters of the z-separation, 9 XZ , and 9 yz varied about 
their respective experimental best-guess values. For each 
FEA run, the values of these three parameters were cho- 
sen at random from Gaussian distributions with means 
equal to the respective best-guess values and standard 
deviations equal to the experimental uncertainties, as 
shown in Tbl. ^ The range of the z-separation was 
considered part of the uncertainty in z; the z- value was 
varied about the middle of this range (26 /xm) and half 
of the range (2 /xm) was added in quadrature to the ex- 
perimental uncertainty. Eight values of A were chosen for 
the FEA models. 



C. Unvaried Inputs to the FEA 

In the FEA model, the mass densities of gold and sili- 
con were taken to be the bulk densities. Though the test 
mass was rotated in the plane of the cantilever, this ro- 
tation was not considered in the FEA model; rotation of 
the test mass would only affect the results at the percent 
level (and would lower a(X)). The test mass was mod- 
eled as a prism (50 x 50 x 30) /xm 3 , with a symmetric 
pattern of 52 blocks on the bottom face to approximate 



16 



the curved section of the test mass. This model underes- 
timates the size of both the prism part and the rounded 
part of the test mass (as compared to what is given in 
Tbl. [IJ; since this underestimate will only increase the 
final a(X) bound, the approximation is acceptable. The 
test mass was taken to be tilted at 0.35 rad from the 
horizontal plane of the cantilever. The model of the test 
mass in the FEA code did not account for the recessed 
face discussed in Sec. Ill Bl Considering the exponential 
dependence of the Yukawa potential, this missing volume 
(on the side face tilted away from the cantilever) would 
only affect the calculated gravitational force at the 2% 
level. This level of error is small compared to other un- 
certainties in a and may be ignored in this analysis. 

The x-position between the drive mass and the test 
mass was taken to be the experimental best-guess value. 
Both the Newtonian and the Yukawa gravitational forces 
do depend slightly on the x-position. However, this de- 
pendence changes with the tilt 9 XZ and there was no 
clear "worst-case" value for x-position to use in the FEA 
model. To reduce computing time, the model computed 
only the force for a range of 200 fim in y-position. The dc 
force for this reduced range was mirrored to make a range 
of 600 /im from which the ac force was derived, a simpli- 
fication that affected the best-fit a at the 10% level. This 
range of y-positions was centered near the best-guess cen- 
ter of the range of data-acquisition and the same FEA 
model was used to fit to each day of data. Due to possi- 
ble tilt of the drive mass, fitting each day of data to the 
same y-range could also incur errors on the order of 10% 
on the best-fit a(A). 

D. Uncertainty in Bimorph Amplitude 

A bimorph amplitude of 125 fim was considered in the 
FEA model. To account for the experimental uncertainty 
of 10 fim in this factor, each Monte Carlo result was also 
fit to the data for bimorph amplitudes in the range of 
95-155 /im (in steps of 5 /an), accounting for three times 
the experimental uncertainty. With the dependence of a 
on the bimorph amplitude thus determined, each best- 
fit value of a for a bimorph amplitude of 125 /im was 
varied 50 times according to a random sampling from 
a Gaussian distribution of bimorph amplitudes, yielding 
a set of 320 x 50 best-fit values of a for each chosen A 
value. Different samplings of bimorph amplitudes (from 
the Gaussian) were used for each day of data. 

E. Uncertainties in Multiplicative Factors 

Several factors in the FEA model were not varied as 
Monte-Carlo parameters. Most of these parameters en- 
tered as multiplicative factors in the determination of the 
best-fit a(X). 

Though the model on which the fit is based (Eq. EJ) 
considers measurements at each y-point, due to thermal 



noise, data will never match FEA predictions at the pre- 
dicted zero minima of the force magnitude. Thus, the 
best-fit a is mostly determined by the quality of the 
fit of the data to the maxima of the prediction curves. 
Roughly, a = (F m — F^)/Fy, where only the maxima 
of the measured force (F m ), the Newtonian force (Fat), 
and the Yukawa force (Fy) are considered. This argu- 
ment holds true in the presence of a small offset force. 
In this experiment, F m 3> FV and thus, the propor- 
tional error in a was determined by a quadrature sum 
of the proportional errors in Fy and F m : (5a/a) 2 ~ 
(5F m / F m ) 2 + (5Fy I Fy) 2 , where 8a is the uncertainty in 
the best-fit a resulting from the uncertainties SF m and 
SFy. This is true for small (< 10%) relative uncertainties 
in Fy and even for large (~ 50%) relative uncertainties 
on F m . The uncertainty in Fn was ignored because the 
contribution of Fn in this equation is small. 

The uncertainties in the measured force (listed in 
Tbl. IIIip were all multiplicative factors, which would scale 
the entire curve of F m (y) up or down. Uncertainties in 
the respective densities of the masses were multiplicative 
factors in Fy. However, all densities were taken to be 
the given bulk densities. Voids and uncertainties in the 
shape of the test mass all entered, to first order, as mul- 
tiplicative factors in Fy . The indeterminate composition 
areas at the bottom of the drive mass were also approxi- 
mated as a very small multiplicative factor on Fy; FEA 
models showed that even decreasing the height of the 
drive mass from 100 /im to 90 fiiii had at most a few per- 
cent effect on the magnitude of Fy for the studied values 
of A. Together, these factors summed to 8.4% relative 
uncertainty on Fy. 

To account for these uncertainties in the multiplicative 
factors of F m and Fy, each of the prediction curves and 
the data curves could be varied by Gaussian distributions 
representing the uncertainty of the respective multiplica- 
tive factors and the best-fits could be sought between the 
new, much larger sets of prediction and data curves. This 
parametric bootstrap was simplified because, as argued 
above, the best-fit a changed with scalings of F m (y) and 
Fy(y) in a predictable way. Thus, the uncertainties in 
these multiplicative factors were counted by varying the 
best-fit a(X) values according to the relative uncertain- 
ties in F m and Fy. 



F. Statistical Uncertainty 

Though a least-squares fit does find a best fit that ac- 
counts for the possibility of fluctuations in the data about 
the "true" values, the statistical uncertainty in the mea- 
surements does lead to an error in the best-fit results. In 
the simple case of a least-squares linear fit, for example, 
this uncertainty in the fit parameters can be easily cal- 
culated [3(J. In this case, the uncertainty in the best-fit 
results cannot be determined analytically. To determine 
the uncertainty in a resulting from the statistical uncer- 
tainty (due to thermal noise) on each measurement, the 
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TABLE IV: Errors on best-fit a 



Parameter 


% Uncertainty 


Volume of test mass 


7.6 


Voids in test mass 


2.5 


Voids in drive mass 


z.o 


Drive mass polished Au/Si boundary 


1.0 


Total Multiplicative error in Fy 


8.4 


Multiplicative error on measured force 


17 


Statistical error 


22 


Total Uncertainty in a 


29 



measurements were artificially varied by the statistical 
uncertainty found in the data. At a given y-point, an off- 
set, drawn at random from a Gaussian distribution with 
a mean of and a standard deviation equal to the stan- 
dard error on Fr, was added to the measured Fr. The 
measured Fj was similarly perturbed. At each y-point, 
the measurements were dithered in this way, adding to 
the measured points their statistical uncertainty. The re- 
sults of a typical FEA simulation were then fit to these 
dithered points. The process was repeated multiple times 
and from this, the uncertainty in a due to the statisti- 
cal uncertainty in the data was found to be 22%. This 
compared well to the standard deviation of the best-fit a 
over many sets of simulated data. 

G. Summation of Errors 

Uncertainties in the most important geometric factors 
(vertical separation and tilt angles) were considered in 
the Monte Carlo run. Uncertainties in the bimorph am- 
plitude were considered by varying the best-fit a(A) re- 
sults. Uncertainties in the multiplicative factors and the 
statistical uncertainty in the best-fit results were consid- 
ered by varying the best-fit results once more. In this 
case, each best-fit result was multiplied by a set of 50 
random samples from a Gaussian of mean 1 and a stan- 
dard deviation equal to the quadrature sum of the uncer- 
tainties due to the multiplicative factors and the statis- 
tical uncertainty on a. These factors are summarized in 
Tbl. |TV| Finally, for each value of A, there were 800 000 
best-fit a results. 



H. Best-Fit Results 

Best- fit yo, 9q, Rq, and Iq in these fits have no im- 
portant physical meaning. However, examination of the 
best-fit results for these other fit parameters was useful in 
evaluation of the analysis method and demonstration of 
the robustness of the fitting procedure. Fig.ll2lshows cor- 
relations among these parameters and the mean square 
error from the best-fit. As expected, the sets of best-fit 




-40 -20 -4 -2 
Best-Fit y Q (urn) Best-Fit R Q (aN) 



FIG. 12: [color online] Correspondence between parameters 
for the fits of the Monte Carlo results to Day 1 data. In 
the upper left, the best-fit a values are shown; from the top 
(larger a) , the groups of best-fit values correspond to values 
of A of 4, 6, 10, 18, 34, 66, 130, and 258 /Ltm. Different symbols 
are used for each value of A, with the same symbols used in 
each plot. Plots show the mean square error (MSE) for the 
fit and the best-fit parameters yo, 6o, Ro, and lo- 



af for each value of A showed an exponential dependence 
on the effective vertical separation (including both z and 
the additional separation resulting from 9 xy and 9 yz ) be- 
tween the masses. 

The mean best-fit yo for each of the two data days in- 
dicated ~ 40 (im difference (for a given CPS reading of 
y-position) between the two days. The best fit 9q var- 
ied by almost tt between the two data days; this is the 
expected change if there were indeed a gravitational sig- 
nal, given the difference in yo. The magnetic experiment 
of Cooldown 02 showed that a shift of this size in the 
CPS reading from one day to the next could be expected. 
Moreover, if the measured force was a non-gravitational 
background (such as electrical noise), an arbitrary shift 
in the best-fit yo between data days could be expected. 

Scaled and shifted by the best-fit parameters, the FEA 
results were almost all within two standard deviations of 
the measured Fr, Fj, force magnitude F, and measured 
phase with respect to the drive. A typical best-fit result, 
in comparison to the measurements, is shown in Figs. 11^1 
andll4l 

For each of the two data days, the mean of the best- 
fit offset force (described by i? and Io) had a magnitude 
close to the mean magnitude of the measured force across 
the day's y-scan. The magnitude of the mean best-fit 
offset for Day 1 was 3.1 x 10~ 18 N, at a phase of 2.4 rad. 
For Day 2, the mean best-fit offset had a magnitude of 
6.3 x 10~ 18 N, at a phase of -2.3 rad. 

The measured force on Day 2 was larger than the mea- 
sured force on Day 1; correspondingly, the typical best-fit 
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FIG. 13: [color online] The measured force as Fr, Fi, mag- 
nitude, and phase for Day 1 as compared to a typical best-fit 
Monte Carlo result, for A = 18 /im. The force magnitude F 

is \ Fr+Fi . Error bars show the statistical error on the 
mean. 



inated from a gravitational coupling between the drive 
and test masses. If the above results were to be inter- 
preted as a signature for a true force, the uncertainties 
would be considered in the opposite sense (with respect 
to the best- guess values) than for this analysis; this would 
yield smaller values of cv(A), a lower bound. Such an ap- 
proach will not be explored in this paper since more data 
with a better signal-to-noise ratio is needed. Future work 
to this end is discussed in Sec. IXII Furthermore, to con- 
firm the existence of a gravitational force, measurements 
would have to be made as a function of the mass or the 
separation between masses in the experiment. 

In the meantime, we will use our current results to 
put new upper bounds on Yukawa-type deviations from 
Newtonian gravity. This is described in the next section. 
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FIG. 14: [color online] The same information as in Fig. 1131 
except for Day 2 data. 

a for Day 2 was 20%~40% higher than the best-ht a for 
Day 1. The mean square errors for the fits on Day 2 
were 20% (at small A) to 40% (at larger A) higher than 
the errors on the same fits for Day 1. 

Though the parameter yo accounted for some of the un- 
certainty in the location of the range over which measure- 
ments were made, the least squares fit did not account 
for uncertainty in the separation between y-points. Un- 
certainty in the n-position between masses was also not 
considered. Furthermore, the possibility of scaling errors 
in F m that changed over the course of the day (such as 
a continual drift of /o away from 3/d over the course of 
the day, changing both the scaling and the phase of any 
measured force) were not included in the analysis. Inclu- 
sion of these effects, all expected to be relatively small, 
has been left for future work. 



I. Interpretation of the Results 

The Monte Carlo results and subsequent varying of the 
best-fit results provide a spread of a for each A consid- 
ered. These results may point to a force with a finite 
magnitude that shows gravity-like features in our exper- 
iment. However, due to the small size of the measured 
force as compared to thermal noise, we are not able to 
determine at this point whether the measured force orig- 



IX. DRAWING A BOUND ON a(A) 

A. Summary of Procedures 

The a(A) bound was derived from a measurement 
of cantilever motion, as described in the previous sec- 
tions. The drive mass was oscillated underneath the 
cantilever bearing the test mass. An interfcrometric 
measurement of cantilever motion was recorded as time- 
series data. The time-series data from the interferome- 
ter were Fourier-transformed, averaged, and compared 
to the drive signal in order to determine how much 
the cantilever was moving at the frequency of inter- 
est (Eq. |5J). From this motion, the force on the can- 
tilever was deduced, as described in Eq. [31 The equi- 
librium position of the drive mass with respect to the 
test mass was varied longitudinally and measurements 
were recorded as a function of this position; room tem- 
perature alignment and capacitive sensors used during 
the low-temperature measurement indicated the relative 
position between masses. The measured force as a func- 
tion of the drive mass position was fit to a model of the 
gravitational force (including a Yukawa force) between 
masses; this model was based on FEA using Eqs. and 
Experimental uncertainties were considered via Monte 
Carlo variation of the inputs to the FEA and variation 
of the results of the fit. The best fit between the mea- 
surements and each FEA result yielded a set of best-fit 
a(X) for a Yukawa potential that could be consistent with 
our data. The final upper-bound results were determined 
from this set of best-fit a(X). 

B. Determination of the Final Results 

The goal is to set an upper-bound on Yukawa-type de- 
viations from Newtonian gravity at the 95% confidence 
level. Though all experimental uncertainties were as- 
sumed to be Gaussian, the resulting distribution of a 
for each A was not; especially in the case of small A, 
where the best-fit a had a highly nonlinear dependence 
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FIG. 15: [color online] An example histogram of the fully- 
varied best-fit q results for A = 18 /im and Day 1 data. The 
Scaled a is a/a and the Relative Frequency is scaled by the 
total number (800,000) of varied best-fit a values. 



TABLE V: Summary of results 



Day 1 



A (pm) MC Mean a a 


Full Set Mean a b 


Upper Bound a c 


4 


3.6 x 10 7 


3.8 x 10 7 


1.5 x 10 s 


6 


6.7 x 10 5 


6.9 x 10 5 


2.0 x 10 6 


10 


2.2 x 10 4 


2.3 x 10 4 


5.1 x 10 4 


18 


2.0 x 10 3 


2.1 x 10 3 


3.7 x 10 3 


34 


5.0 x 10 2 


5.1 x 10 2 


8.4 x 10 2 


66 


2.7 x 10 2 


2.8 x 10 2 


4.4 x 10 2 


130 


2.2 x 10 2 


2.3 x 10 2 


3.6 x 10 2 


258 


2.1 x 10 2 


2.1 x 10 2 


3.3 x 10 2 


Day 2 


A (/im) 


MC Mean a 


Full Set Mean a 


Upper Bound a 


4 


5.1 x 10 7 


5.2 x 10 7 


2.0 x 10 s 


6 


8.9 x 10 s 


9.1 x 10 5 


2.8 x 10 6 


10 


2.9 x 10 4 


3.0 x 10 4 


6.8 x 10 4 


18 


2.5 x 10 3 


2.6 x 10 3 


4.8 x 10 3 


34 


6.2 x 10 2 


6.4 x 10 2 


1.1 x 10 3 


66 


3.4 x 10 2 


3.5 x 10 2 


5.7 x 10 2 


130 


2.7 x 10 2 


2.8 x 10 2 


4.6 x 10 2 


258 


2.5 x 10 2 


2.6 x 10 2 


4.2 x 10 2 



"This is the mean of the best-fit a results from fitting the set of 
Monte Carlo-varied FEA prediction curves to data. 

''This includes the variation of Monte Carlo results to account 
for uncertainty in bimorph motion and the statistical uncertainty 
in a combined with experimental uncertainty in the multiplicative 
factors. 

c This is the 95% confidence level result. 



on the vertical separation and tilts between masses, the 
distributions were very asymmetric. No assumption was 
made of the analytic form of the distribution; the 95th 
percentile a from the distribution provided the result at 
the desired (one-sided) confidence level, as diagrammed 
in Fig. E| The 95th percentile cv(A) after the variations 
of the Monte Carlo results was found to be less than 35% 
more than the 95th percentile a(X) from the Monte Carlo 
variation alone. 

Though this analysis assumed a > 0, because the mea- 
sured force was ~ 200 times greater than the Newtonian 
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FIG. 16: [color online] The parameter space of cv-A of Yukawa- 
type deviations from Newtonian gravity. Results from this 
paper (Stanford 2) are shown in comparison to previous re- 
sults from this experiment js| (Stanford 1). Also included 
are theoretical predictions (dashed line) and an experimental 
result (solid line) as shown Fig. 



gravitational force in this case, the results for a < 
would be statistically indistinguishable (when compar- 
ing |a|) from the results for a > 0. Thus, our results 
may be considered as a bound on \a\. 

Results for Day 1 were different from results for Day 2. 
As discussed above, there was indication of a larger non- 
gravitational background on Day 2 as compared to Day 1. 
This could be the cause of the slightly higher best-fit 
a(A) on Day 2. If the difference between the mean best- 
fit a(A) and the 95th percentile a(X) is considered to 
be two standard deviations, then the difference between 
Day 1 and Day 2 results is at an acceptable level. The 
Day 1 results are plotted in Fig. ^| and taken as our 
95% confidence-level bound. Results from Day 1 were 
chosen because the error on the fit and diagnostics from 
the data (such as signal levels at fa) suggested that Day 1 
measurements suffered from a smaller non-gravitational 
background. 

These results are properly interpreted as a bound on, 
not a discovery of, deviations from Newtonian gravity. 
The choice of the Yukawa parameterization is appropri- 
ate though not all-inclusive; results could be similarly 
analyzed using other parameterizations (such as a power- 
law potential) for deviations from Newtonian gravity. We 
are 95% confident that no Yukawa-type gravitational po- 
tential exists with cv(A) above the bound reported for the 
Day 1 results in Tbl. [7| 
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X. CONCLUSIONS 

In this paper, we have described our experimental ap- 
paratus and the latest data it has produced. While this 
experiment was based on the same principles as the one 
that yielded our first results Q, it incorporated many 
substantial improvements in the design, technique, and 
data analysis. The results presented here represent al- 
most an order of magnitude improvement over previous 
results at A ~ 20 /im |9j, yielding the most stringent ex- 
perimental constraints on Yukawa-type deviations from 
Newtonian gravity at length scales of 6-20 jam. This 
new bound provides constraints on predictions of mod- 
uli and gauge bosons. These new constraints do not rule 
out string theory or supersymmetry or the possibility of 
large extra dimensions. However, the a(X) bound limits 
what kind of particles could be included in any of these 
theories of physics beyond the standard model. 

XI. FUTURE PROSPECTS 

These results are a quantitative and qualitative im- 
provement over previous results. Moreover, our appara- 
tus and analysis are more complete and robust, provid- 
ing a strong foundation on which to build future exper- 



iments to test gravity even further. Reducing electrical 
and mechanical coupling between the bimorph and the 
interferometer and using a test mass without a curved 
face (to improve fiber-mass alignment, increasing fringe 
height) could increase the signal-to-noise by an order of 
magnitude. It is unlikely that the thermal noise limit 
in this experiment will be reduced by more than a fac- 
tor of two. However, measuring with a test mass not 
tilted with respect to the cantilever would increase the 
expected gravitational force and could easily yield a fac- 
tor of two improvement on the a(X) bound. A new can- 
tilever design may provide a magnetic calibration, which 
will greatly reduce uncertainty in y-position and allow 
for a more precise fitting of results to FEA predictions. 
We are developing a similar apparatus with a circular 
drive geometry, which we expect to be at least an order 
of magnitude more sensitive than the current one. 
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